Role of Polaron Pair Diffusion and Surface Losses in Organic Semiconductor Devices 
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By applying Monte Carlo simulations we found that the extraction of bound polaron pairs (PP) at the elec- 
trodes is an important loss factor limiting the efficiency of organic optoelectronic and photovoltaic devices. 
Based upon this finding, we developed a unified analytic model consisting of exact Onsager theory, describing 
the dissociation of PP in organic donor-acceptor heterojunctions, the Sokel-Hughes model for the extraction of 
free polarons at the electrodes, as well as of PP diffusion leading to the aforementioned loss mechanism, which 
was not considered previously. Our approach allows to describe the simulation details on a macroscopic scale 
and to gain fundamental insights, which is important in view of developing an optimized photovoltaic device 
configuration. 
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In solar cells based on polymer-fullerene bulk hetero- 
junctions (BHJ) primary molecular excitations, also called 
excitons, are generated by light absorption. The photo- 
generated excitons diffuse by emission and re-absorption to 
donor-acceptor interfaces.^ There, the excitations may un- 
dergo charge transfer from the excited donor molecules to the 
nearby acceptors. 23 The resulting charge transfer states are 
also called polaron pairs (PP) and, unlike bipolarons, consist 
of oppositely charged constituents. They may dissociate into 
free polarons which can be extracted from the device to gen- 
erate a net photocurrentFl In organic solar cells the net pho- 
tocurrent is crucially determined by the PP dissociation and 
extraction yields, currently leading to overall power conver- 
sion efficiencies of up to 8%P 

It is commonly implied that PP either recombine or disso- 
ciate at the place of their creation, as their diffusion is not 
accounted for in theoretical models. 6 8 However, it is known 
from experiments that photogenerated charges have very high 
initial mobilities, 9 and there are indications that photogener- 
ated charges can move to different molecules or conjugation 
segments before they recombine or dissociate. 10 In this Letter 
the diffusion of Coulomb bound PP resulting in an important 
loss mechanism in organic BHJ solar cells is considered. We 
perform Monte Carlo simulations of PP recombination and 
photocurrent extraction in BHJ, focussing on PP diffusion and 
the resulting losses at the electrodes. We find that without 
blocking layers, over 40% of the generated PP can get lost at 
the electrodes and that, for high electric fields, this loss mech- 
anism becomes dominant. Such high surface loss cannot be 
explained by the theories commonly used. Our simulations 
show that PP diffusion is responsible, as it can increase ex- 
traction of both (bound) charges at the same electrode. Upon 
this finding we present a unified analytic model to calculate 
the photocurrent yield, accounting for PP diffusion and losses 
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at the surface. 

In our Monte Carlo simulations we describe the organic 
blend of P3HT (poly(3-hexyl thiophene)) and PCBM ([6,6]- 
phenyl-C6i butyric acid methyl ester), a widely studied 
representative of polymer-fullerene semiconductors, on a 
mesoscopic scale. 11 Therefore, electron donor and acceptor 
molecules were distributed in a volume ratio of 1 : 1 within 
a 100 x 25 x 25 cubic lattice with constant spacing of 1 nm, 
sandwiched between electrodes on the long lattice side and 
periodic boundary conditions for the other directions. Each 
acceptor molecule was assigned to a single cubic lattice site, 
in accordance with the spherical shape of fullerenes. How- 
ever, as conjugated polymers consist of several monomer units 
and are known to have an effective conjugation length (CL), a 
donor molecule was spanned over multiple lattice points.^ 

Charge tra nsport was described by the Miller- 
AbrahamfjEE] rate equation, also accounting for charge 
carrier Coulomb interaction as well as mirror charge effects. 
Miller-Abrahams theory was used as it neglects all polaronic 
effects and thus qualitatively better accounts for PP excess 
energy, which is initially (over-)compensating molecular 
reorganization energies in e.g. classical Marcus theory. 
In order to include energetic disorder, the energy levels, 
corresponding to the respective molecular orbitals, were 
taken from Gaussian distribution functions according to our 
experimental findings, with standard deviation Ga — 60 meV 
for acceptor molecules and Go = 75 meV for donors.^ The 
PP lifetime was set to x e ff = 10~ 7 s! 11 * 14 * 15 ! At simulation 
initialization, one polaron pair was set at a randomly chosen 
donor-acceptor interface, which corresponds to a PP density 
of npp = 1.6 • 10 16 cm~ 3 . To gain statistically reliable results, 
about 2.5 ■ 10 s simulations were performed for each parameter 
set. 

A crucial step during the photocurrent generation is the PP 
dissociation yield, /?,■„,. We point out that p,„, corresponds 
to internal polaron pair dissociation, where surface effects 
and losses are neglected. We considered it in our former 
workEl in which charge carrier delocalisation, implying high 
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Figure 1: (Color Online) Probability of PP surface losses vs. 
electric field found by Monte Carlo simulation of 1 : 1 donor- 
acceptor blends at 300 K with % e tt = 100 ns (device length 
/ = 100 nm, dielectric constant £ = 3.0). (a) Absolute prob- 
ability for PP surface loss pi oss is significant at low fields and 
for long donor conjugation lengths (CL). (b) pi oss relative to 
the sum of all losses — in the bulk (p rec ) and at the surface — 
becomes the dominant loss mechanism for increasing CL and 
increasing field. 



local mobilities, was found to be a key parameter for achiev- 
ing high pi, lt . In this Letter, we consider losses at the elec- 
trodes, and therefore define the charge extraction yield, p ext , 
which is a measure for the photocurrent. We show that by in- 
creasing CL, the probability for losing PP at the electrodes, 
Pioss = Pint — Pext, increases (Fig. [lj). Even more, the frac- 
tion of PP losses at the surface pi oss becomes dominant for 
CL=10 over all electric fields, as compared to the sum of all 
bulk {p rec = 1 — p^t) and surface PP loss mechanism. 

Investigating this significant loss mechanism in detail 
showed that the PPs do not necessarily dissociate or recom- 
bine at the place of their creation. Due to their strong mutual 
Coulomb binding, the generated polaron pairs diffuse as neu- 
tral quasiparticles within their lifetimes along the distributed 
donor-acceptor interfaces. To illustrate this effect, the diffu- 
sion of external field stabilized PPs in direction of that ex- 
ternal field of 3 ■ 10 5 V/m is shown in Fig. [2] exemplarily for 
CL = 6. The graph shows the fraction of PP vs. their diffu- 
sion length before recombination occurs. Only a tenth of the 
PP does not diffuse, as indicated by the peak of the probability 
distribution. For longer CL, the PP diffusion length increases, 
whereas higher external fields reduce the PP bulk recombi- 
nation in general (not shown). Thus, the surface loss of PP 
increases at the same time as bulk recombination p rec reduces 
with increasing CL.^ 

As the diffusion lengths distribution was averaged over 
more than 10 5 simulations, the discrete random-walk dif- 
fusion process of the simulation can be approximated by a 
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Figure 2: (Color Online) Probability distribution po pp of polaron 
pair diffusion lengths from Monte Carlo simulation (circles) of 
a 1 : 1 donor-acceptor blend at 300 K with % e ff = 100 ns (de- 
vice length / = 100 nm, dielectric constant £ = 3.0). The dashed 
lines indicate the fitted contributions of the two identified diffu- 
sion processes, with diffusion lengths of 2 resp. 9 nm, the solid 
line representing the sum of both. 



continuous -time stochastic Wiener process, also known as 
Brownian motion. For PP, decaying exponentially with l e ff, 
the probability to diffuse for the distance x is given as 



Pd pp (x)= / exp - 
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2 IDppt 



^ dt . 
(1) 

Thereby Dpp is the PP diffusion constant, which relates to 
the characteristic diffusion length lj of the process as Id — 

y/Dpp-X e ff . 

We find that the simulated distribution can be expressed 
well as the combination of two diffusion processes of differ- 
ent Dpp, as shown in Fig. [2] For the given parameters, 40% 
of the PP hardly diffuse and recombine within a distance of 
l c i = 2 nm. The other 60% diffuse for l c ( — 9 nm, in total cor- 
responding to 12% surface loss due to PP diffusion (device 
length 100 nm). For thinner devices, the relative losses are 
even more significant (not shown). We observe that the two 
diffusion processes found are related to the energetic relax- 
ation of charge carriers in the initially equally occupied den- 
sity of acceptor or donor states. 

To study these high surface losses found in our simulations 
in more detail, and to identify the essential contributing pro- 
cesses, we include the PP diffusion into an analytic model. 

Our approach is shown schematically in Fig. [3] and is as 
follows: Starting from the created PP, first their change in 
number is described by competing rates (Fig. [3] left triangle), 
similar to Onsager-Braun 6 7 theory. However, by not restrict- 
ing our approach to time independent rate constants, we can 
include a time dependent surface loss process for the photo- 
generated PP. In contrast, Onsager-Braun or exact Onsaget^Sl 
extensions are entirely done on infinite space, thus c ompl etely 
neglecting surface losses due to space confinement!^^ Sec- 
ond, the final extraction yield of the individual charge carriers 
created by PP dissociation is calculated with a drift and dif- 
fusion model for the individual charge carriers (Fig. [3] right 
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Figure 3: (Color Online) Combination of dissociation theory (left 
triangle) and drift and diffusion model (right triangle). PP dis- 
sociation probability Pdiss results from competitive rates of the 
individual PP loss processes Kj and leads to net current with 
Psh- Direct surface loss of PP (Ki oss ) has not been addressed 
before. 



triangle)PEl 

In our general, time dependent approach, the probability 
Pi of a PP to decay by a certain process i — such as recombi- 
nation, dissociation or surface loss — is the integral over the 
relative probability k,(t) of that individual decay process and 
the probability n(t) that the PP still exists at the time t after its 
creation, 



Pi = / ki{t)n{t)dt 
Jo 



(2) 



The time dependence of kj(t) is relative to the creation of the 
PP and not due to external changes on the whole system, thus 
equal for each PP. Time independent ki are also often called 
(decay) rates and — assuming non-interacting PP initially cre- 
ated at t = Os — n(t) is called the number of (existing) PP. 

The competing, time dependent decay probabilities ki(t) re- 
sult in the general differential equation for n(t ), 



n{t) = -'£ki(t)-n(t) 



(3) 



where the formal solution for the initial value of n(0) = 1 can 
be written as 



n(t)=Y[Ki(t) with Ki{t) = 



(4) 



We note that for time independent kj the integral over n(t) 
in Eq. Q results in the inverse sum of ki and so yields the 
common equation of competing rates, as e.g. known from 
Onsager-Braun theory. 

In order to calculate the probabilities of the PP losses p, 
according to Fig. [3](left triangle), ki(t) or equivalent K{(t) for 
recombination, surface loss and dissociation process have to 
be determined. 

Commonly, recombination processes are described accord- 
ing to exponential decay with a recombination rate k rec , which 
is associated with the inverse PP lifetime x by k rec = 1 /x. 

In order to calculate the PP dissociation rate &</to> we 
start from an exact solution of Onsager theory^ to get an 
expression for the PP dissociation probability which con- 
tains recombination and dissociation losses, but neglects any 
surface effects. In detail, the PP dissociation probability 



tydiss(krec,E ,rpp,p) accounts for Langevin type recombina- 
tion, finite recombination rates and finite recombination dis- 
tances and mainly depends on recombination rate k rec , exter- 
nal field E, initial PP radius rpp and charge carrier mobility 
^uPln a second step, the dissociation rate k c n ss is obtained out 
of the dissociation probability ffdiss by inverting its equivalent 
expression as competing rates of k rec and kdi SS > thus 



kdiss 



tydiss ■ 



1 - tydiss 



(5) 



Alternatively, kdiss could be approximated directly with 

Onsager-Braun!- 7 -' 

In the newly considered surface loss process a PP is lost 
when it diffused to one of the electrodes. As the PP are cre- 
ated uniformly within the device, surface loss is the dominant 
loss mechanism at early times because PP generated close to 
the electrodes most likely get lost. However, at longer times 
surface loss rapidly becomes less significant (not shown), as 
the PP loss probability changes during the lifetime of the PP. 
In order to describe the time dependent surface loss process 
due to PP diffusion with the diffusion coefficient Dpp, the 
same continuous -time stochastic Wiener process as in Eq. ([T]) 
is used, 



1 H ! H' 1 / 1 1 -\2 
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dxdx. 



(6) 



The real device length has to be taken relative to the average 
hopping distance of the charge carriers, because the Wiener 
process is the continuous limit of a discrete random walk pro- 
cess. We express this by a virtual device shortening to an ef- 
fective device length l e tt. Thus, the importance of the surface 
loss process is proportional to the ratio of diffusion length l c / 
to the effective device length l e ff. 

Now, p c iiss and pi„ r = 1 — p rec can be calculated by inte- 
grating Eq. (|2| with k rec , kdiss an d K[ oss (t) from above, which 
is the last step in calculating the net generation of free charge 
carriers (Fig. [5] left triangle). We note that for ki oss (t) = 0, i.e., 
Ki ss(t) = 1, Eq. (|2| simplifies to exact Onsager or Onsager- 
Braun theory, depending on which theory was used to describe 

kdiss- 

Once a PP is dissociated into free polarons, these free 
charge carriers have to be extracted over the electrodes to con- 
tribute to the net photocurrent (Fig. [3] right triangle). Unlike 
the pure PP diffusion of neutral quasiparticles, the individual 
free charges are affected by the external field. 

A model for non-interacting charge carriers describing 
charge extraction in a device due to drift and thermal diffu- 
sion is that from Sokel and Hughes. 17 Assuming a uniform 
charge carrier generation without recombination losses within 
the device, the probability p co u of collecting charges with an 
external electric field E on the corresponding electrodes of a 
device of real lenght / is 
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exp(eEl/k B T) + l 2k B T 
exp(e£7 /k B T) — 1 eEl 



(7) 
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Figure 4: (Color Online) Probability for net photocurrent p ext and 
PP dissociation p ln( from Monte Carlo simulation (markers; 
1 : 1 donor-acceptor blend, CL= 6, X e tf =100 ns, 300 K, / = 
100 nm and e = 3.0) and analytic calculations (fit; Dpp = 3.8 • 
10~ 10 m 2 /s, l eff = 22 nm, r PP = 3.0 nm, p = 4.3 • 10~ 8 m 2 /Vs , 
X e ff = 100 ns, 300 K, / = 100 nm and e = 3.0). For comparison, 
calculations without PP diffusion (Dpp = m 2 /s) and without 
virtual device shortening (leff = 100 nm) are shown. 



The probability of losing individual, free charges on the wrong 
electrodes is simply 1 — p co li (Fig. [3] right triangle). 

The total probability of net photocurrent generation is fi- 
nally given as the product of net generation and collection 
probabilities as p ext = p C(M ■ p diss . 

The application of this new theoretical approach to explain 
our simulation results is shown in Fig. [4] Wherever possi- 
ble, the parameters for the calculations were taken equivalent 
to those in our simulation, in particular X = l e ff = 10~ 7 s, 
T = 300 K and I = 100 nm. Dpp was also set to a value de- 
termined in our simulation, i.e., to Dpp = 3.8 ■ 10~ 10 m 2 /s 
as average of the two diffusion processes (Fig. |2j. rpp, fj 
and l e ff were taken as fitting parameters. Best agreement 
was found for rpp — 3.0 nm and /j = 4.4 ■ 10~ 8 m 2 /Vs, val- 
ues that are in good agreement with our recent publicatiorP=i 
For l e ff = 22 nm the device length of 100 nm is shortened 
virtually by a factor of 4.5 for CL= 6. We find perfect agree- 
ment between our simulation results and the extended theory 
(Fig.g). 

For comparison, calculations were also done for Dpp = 



m 2 /s and for l e ff = 100 nm. Clearly, the simulation results 
cannot be reproduced with these parameter sets. 

In conclusion, by using mesoscopic Monte Carlo simula- 
tions of polaron pair dissociation and charge extraction in 
polymer-fullerene blends, we found that the higher local mo- 
bility directly increases the polaron pair diffusion length, and 
delocalisation leads to less hops being needed to overcome a 
given device length. Thus, the probability of neutral Coulomb 
bound polaron pairs to accidentally diffuse to an electrode 
within their given lifetime grows with increasing delocalisa- 
tion length. In the framework of our Monte Carlo simulation, 
these effects lead to a high loss of polaron pairs at the surface, 
thereby reducing the photocurrent. 

In order to describe our findings analytically, we presented 
a unified model consisting of exact Onsager theory for po- 
laron pair dissociation, the Sokel-Hughes model for charge 
extraction, and our extension to account for losses by polaron 
pair diffusion. Applied to organic semiconductor devices, this 
extended model leads to a deeper insight into the relation of 
local morphology and the photocurrent. This has important 
consequences for the optimization of organic optoelectronic 
and photovoltaic devices, highlighting the need to reduce sur- 
face losses by adjusting the device configuration in view of 
blocking layers or selective electrodes. Indeed, in bulk hetero- 
junction solar cells with conjugated polymer donors inhibiting 
high conjugation lengths, blocking layers not only suppress 
the extraction of free polarons at the wrong electrodes. They 
also minimize the diffusion of Coulomb bound polaron pairs 
to the electrodes, thus increasing their dissociation and extrac- 
tion probability. 
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